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Abstract 



We discuss some properties of the exact four-spinon dynamical correlation 
function S4 in the antiferromagnetic spin 1/2 AAA-model, the expression 
of which we derived recently. We show that the region in which it is not 
identically zero is different from and larger than the spin-wave continuum. We 
describe its behavior as a function of the neutron momentum transfer k for 
fixed values of the neutron energy to and compare it to the one corresponding 
to the exact two-spinon dynamical correlation function S% . We show that the 
overall shapes are quite similar, even though the expression of S4 is much 
more involved than that of S^- We finish with concluding remarks. 
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I. INTRODUCTION 



Quantum spin chains have been the subject of intensive study during the past seven 
decades |]J . Experimentally, their properties are investigated via inelastic neutron scattering 
on (anti)ferromagnetic quasi one- dimensional compounds 0. One important quantity that 
holds much of the information related to such compounds is the dynamical correlation 
function (DCF) S of two local spin operators. Indeed, the neutron scattering intensity 
is directly proportional to it, see for example @. 

An important feature of these chains is that some of them, like the Heisenberg model, 
are amenable to exact theoretical treatment while still describing nontrivial interactions [|TJ , 
see also 0. This is because they incorporate in them a rich mathematical structure: the 
quantum affine algebra. The early work on such models consisted in determining exactly 
their static thermodynamic properties, see |IJ, whereas, more recently, a fuller exploitation 
of the quantum symmetry using bosonization techniques has allowed for a more systematic 
description of their dynamical properties 0] . A systematic account of this work is given in 

B- 

However, one aspect of these new techniques is that the exact correlation functions are 
usually obtained in the form of quite complicated contour integrals and this renders their 
manipulation somewhat cumbersome. But on the other hand, before this exact treatment, 
the approach to these dynamical quantities was only approximate. Indeed, if one considers 
for instance the DCF in the antiferromagnetic Heisenberg model, the focus has for a long 
time been only on what we now know to be its two-spinon contribution S2. First there has 
been the Anderson (semi-classical) spin- wave theory ||, an approach based on an expansion 
in powers of 1/s, where s is the spin of the system and hence, is exact only in the classical 
limit s = 00. It can describe with some satisfaction compounds with higher spins [[F], but 
fails in the quantum limit s = ~. Then there has been the so-called M filler ansatz ||, 
which gives an approximate expression for S2 that can account for many aspects of the 
phenomenology for s = \ compounds. Only very recently could we get for this system a 



final exact expression for 5*2 0, and that gave a better account of the data [10]. 

Now for the spin-| antiferromagnetic Heisenberg model, there is a need to go beyond 
the two-spinon contribution. The need is two-fold. First, it is important to see if one is 
able to get useful information from these complicated and compact expressions we alluded 
to above. Second and more important perhaps is the fact that, though the exact two-spinon 
contribution accounts for much of the phenomenology as we said, about 70%, in a sense 
that will become clear later in section 3, it still doesn't account for all of it. This point is 
demonstrated in particular in [IT|. 



The natural step forward is to look into the exact four-spinon contribution. To the best 
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of our knowledge, refs |T2[ and [|13j constitute the first direct attempt in this direction. A 
general expression for the n-spinon contribution to the DCF in the anisotropic Heisenberg 
model is given in [O] and (a still compact one) for the isotropic limit in [12]. In ref [13 



we specialize to the four-spinon case and give a discussion of some of its properties in the 
isotropic and Ising limits. In particular, we show that in the isotropic limit, the one of 
interest here, the exact four-spinon contribution is safe of any potential divergences. 

In this work, we further the description of the four-spinon contribution S4. We determine 
the region in which S4 is not identically zero and show that it is larger than that of S%. Also, 
we show that the behavior of S4 as a function of the neutron momentum transfer k is similar 
in its overall shape to that of the corresponding S2. This result is to be contrasted with 



the fact that the expression of S4 is a lot more involved than that of S2, see (^3|) and (|13D 
respectively. 

This paper is organized as follows. In the next section, we briefly discuss the antiferro- 
magnetic spin 1/2 Heisenberg model. We describe the spinon Hilbert-space structure and 
define the dynamical correlation function. In section 3, we succinctly review the properties 
of the two-spinon contribution. We give the exact expression of S% and that of the Miiller 
ansatz, and briefly compare their main features. In section 4, we give the exact expression 
of 64 and discuss in detail its properties we mentioned above. The last section comprises 
concluding remarks. 



II. THE EXACT DCF IN THE 1II-MODEL 

The antiferromagnetic s = \ AAA-Heisenberg model is defined as the isotropic limit of 
the X X Z- anisotropic Heisenberg Hamiltonian: 

1 00 

H xxz = "2 £ «< + i + <<+x + , (1) 

n=— 00 

where A = (q + q~ 1 )/2 is the anisotropy parameter. The isotropic antiferromagnetic limit 
is obtained via the limit q — > — 1~ or equivalently A — > — 1~. Here a^' y,z are the usual Pauli 
matrices acting at the site n. The exact diagonalization directly in the thermodynamic limit 
of the Hamiltonian in (P is possible using the U q (sl(2)) quantum group symmetry present 
in the model f|. The resulting Hilbert space T consists of n-spinon energy eigenstates 
|£i, ■■.,£ n )ei,...,e n ;i built on the two vacuum states |0)j, i = 0, 1 such that: 

Hxxz\£l, Cn)ei,...,e n ;j = £e(£j)|£i, £n)ei,...,e„; i , (2) 

where e(£j) is the energy of spinon j and £j is a spectral parameter living on the unit circle. 
In the above equation, Ej — ±1 and the index i refers to the boundary condition on the spin 
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chain, see |J. The translation operator T which shifts the spin chain by one site acts on the 
spinon eigenstates in the following manner: 

n 

T\£i, •••,6n)e 1 ,...,e„;i = t(&) Cn) ei ,...,e wi l-i , (3) 

i=l 

where t(£j) = e~ tPj and pj is the lattice momentum of spinon j. The expressions of the spinon 
energy and lattice momentum in terms of the spectral parameter are quite cumbersome in 
the anisotropic case, see but simplify considerably in the isotropic limit, see eq ([12]) 

below. The completeness relation in JF reads: 

I=EE E ^ f l6>-^n)e 1) ... ) e n ; ii;ei ,. ie „(6,...,en| • (4) 

i=0,l n>0 {e i= ±l}j =lin ' J j=l ^ 

The two-point DCF is the Fourier transform of the zero-temperature vacuum-to- vacuum 
two-point function, i.e., it is defined by: 

/oo 
dtJ2e i{Mm) ,(01^(^^(0)10),, (5) 
■°° mez 

where u; and k are the neutron energy and momentum transfer respectively and a ± denotes 
(a x ± zo*)/2. The DCF is such that: 

S 1 ' + -(uj, k) = S*'+-(w, -k) = S*' + -{u, k + 2tt) , (6) 

the two relations that express the reflection and periodicity symmetries on the linear chain. 
Inserting the completeness relation and using the Heisenberg relation: 

a x J> z (t) = eMiHxxz t) T~ m a^ z {Q) T m ex P (~tH xxz t) , (7) 

we can write the transverse DCF as the sum of n-spinon contributions: 

fi*+-( W| *) = £ Sil + ~(u;,k), (8) 

n even 

where the n-spinon DCF S n is given by: 



meZ ei,...,e„ " j=l ^ J \ j=l 

relation in which X 1 denotes the form factor: 

K,...,e n (^-^n)= M°a (0)16,-, (10) 
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Note that each S n must satisfy relations @. 

The exact expression of this form factor has been determined in ||. Using this form 
factor, we can give an exact expression for the n-spinon DCF in the anisotropic case, see [|13|| , 



and determine exactly its isotropic limit, see ||12|| . This limit is obtained via the replacement 

IB 



£ = ie- Ae "; q = -e~ £ , e - 0+ , (11) 

where p is the spectral parameter suited for this limit. The expressions of the energy e and 
momentum p in terms of p then read: 

7T 

e (p) = tt^ — r = —Ti" sin p ; cot p = sinh(27rp) ; — n < p < . (12) 

cosh(27rp) 

The transverse two-spinon DCF 5*2 does not involve a contour integration and has been 
given in || . The four-spinon one S4 involves only one contour integration and its expression 
is given in [O. In the next section, we review the properties of S 2 and in the one that 
follows it, we discuss those of S4. 



III. THE EXACT TWO-SPINON DCF 

The exact expression of the two-spinon DCF of the spin | XXX-model is given in || 
and reads: 

1 e" /(p) 

S2 (cu,fc-7r) = - =Q(u -lj 2 i)Q(uj2u-u) , (13) 



4 V^ 



where is the Heaviside step function and the function I(p) is given by: 

dt cosh(2t) cos(4pt) - 1 t 



J ^ /o t sinh(2t) cosh(t) 6 " 

^2u(i) is the upper (lower) bound of the two-spinon excitation energies called the des 
Cloizeaux and Pearson P|| upper (lower) bound or limit. They read: 

k 

uj 2u = 27rsin-; iv 2 i — tt | sin k\ . (15) 

The quantity p = pi — p 2 and is related to uj and k by the relation: 







LU 2 - 





coshnp=J-j^-f, (16) 
a relation obtained using eq (|T2|) and the energy-momentum conservation laws: 
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u = ei + e 2 ; k = -p 1 -p 2 - (17) 

Note that the explicit expression of S 2 given in eq (|i~3"|) satisfies the reflection and periodicity 
symmetries expressed in (|5|). 



The properties of 5*2 have been discussed in ||10|| . There, a comparison with the Miiller 
ansatz |J is given. This latter was derived from the properties of some solutions to the Bethe- 
ansatz equations, from numerical calculations on finite spin chains and from an analysis of 
phenomenological results. It reads: 

where A is a constant determined in such a way to fit best the phenomenology ||. There 
are two main differences between the exact expression fllTf ) and the approximate one (|18|). 
10| . First, the two-spinon threshold at io 2 \ is more singular in ([13|) than in (|18|) . Second, at 



the upper two-spinon boundary uj 2u , S 2 vanishes smoothly whereas has a sharp cut-off. 
But if one defines the frequency moments of the DCF: 

/+00 
duuj n S{uj,k) , (19) 
-00 

one shows that as k — ► 0, the moment of S 2 vanishes as: 

Ki 2 \k)r.^\k), (20) 

and the same holds for the Miiller ansatz . 

Actually, the frequency moments (|1^) are particular cases of a set of sum rules the DCF 
is known to satisfy exactly. For exemple, we know that the first moment is exactly equal to 

0: 



4 In 2 — 1 

Ki(k) = (1-cosJfe). (21) 

It turns out that the same frequency moment for S 2 is such that [1C]: 

- ™% . (22) 
This means that, according to this sum rule, S 2 is way off the total DCF S by roughly 



30%. In fact, other exact sum rules confirm this trend, see [10]. Of course, those remaining 
30% are "filled", so to speak, by the n > 2-spinon contributions. The natural question that 
comes to mind is: how much S4 takes up from these 30%? Obviously, one has first to study 
the behavior of this contribution before trying to answer this question, and this we do in 
the remainder of this paper. 
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IV. THE EXACT FOUR-SPINON DCF 



The analytic expression of the four-spinon DCF has been worked out in [IS] and for 
< k < it it reads: 



/0 rO 
dp 3 / dp 4 F(pi, ...,p 4 ) , 

where C 4 is a numerical constant and the integrand F is given by: 

pi \ = f(Pu-,P*) E?=il&(Pi»-»P4)| s 

The different quantities involved in this expression are: 

W = oj + Ti (sin p 3 + sin p 4 ) ; 

W u = 2vr sin y ; 

K = k + p 3 + p 4 ; 
cot pj = sinh(27T/3j) , 



-7T < Pj < . 



The function / is given by: 



/(pi, ...,p 4 ) = exp 

where pjji = pj — py and the function gi reads 

4 

^=(-l) £+1 (27T) 4 ^COsh (27TP,) 



i<i<i'<4 



X 



~ n^(^ - |g(f - g + igg) A r(m - I + zp„ 



m=eo— 



IL^ 71 " ^mh^ 



J- r(m + 1 + ipji 



(23) 



(24) 



(25) 



(26) 



(27) 



where 9 is the Heaviside step function. In Q2~3|), the sum E( Pl P2 ) i s over the two P arr s (pi,P2) 
and (p2,Pi) solutions of the energy- momentum conservation laws: 



W = — 7r(sin p\ + sin p 2 ) 



They read: 



. K W 

\Pl,P2) = 7T + arccos 7, • K I ' 

2 \ 2tt sin - ' 



K = -pi - p 2 



K I W 

arccos 

2 I 2vr sin f 



(28) 



(29) 
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Note that the solution in (p9|) is allowed as long as Wi < W < W u where W u is given in fl25|) 
and: 

Wi = tt| sin K\ . (30) 

We henceforth put ourselves in the interval < k < n. To get the behavior of S4 for the 
values of k outside this interval, one uses the symmetry relations @ given in section 2. In 
the work |I3|, we have discussed the behavior of the function F given in We have shown 



that the series gi is convergent. We have also shown that in the region where two p^'s or 
more get equal, the function g% is finite. The function / going to zero in these same regions 
| |10| |, this means the integrand F of S4 has a nice regular behavior there. Furthermore, we 
have shown that F is exponentially convergent when one of the pi's goes to infinity, which 
means the two integrals over p 3 and P4 in (^) do not yield infinities. All these analytic 
results pave the way to "safe" numerical manipulations. 

The first thing we wish to discuss in this work is the extent of the region in the (k, uS)- 
plane in which S4 is not identically zero, and compare it to that of S2. Remember that from 
(|T3|), 5*2 is zero identically outside the spin- wave continuum uj 2 i{k) < uj < uj 2u (k), where 
^>2i,u{k) are given in (|l^). From the condition W\ < W < W u discussed after eq (0), one 
infers that in order for 5*4 to be nonzero identically, one has to have uj±i < uj < uj^ u , where: 

^4i{k) = 37rsin(/c/3) for < k < tt/2 ; 

^>Ai{k) = 37rsin(/c/3 + 27r/3) for n/2 < k < n ; 

u iu (k) = 4ttcos(A;/4) for < k < vr . (31) 

All these branches are plotted in fig. 1. 

The first thing we immediately see is that the ^-region, i.e., the region in which S4 
is not identically zero, is not confined to the spin- wave continuum delimited by the dCP 
branches O02i, u given in ([15]). This means that, a fortiori, the full S is also not confined to 
the SVregion. This fact is confirmed by early finite chain numerical calculations || and 
the phenomenology 0. Furthermore, fig. 1 shows that for < k/n < 1/2, the 5*4-region 
is entirely beyond the ^-region. This means that for this interval, we may expect S2 to 
be dominant in 5* within the spin- wave continuum. However, for 1/2 < k/n < 1, there is 
overlap between the two regions such that the ^-region is more or less within the ^-region. 
We may therefore expect here the contribution of S4 to play a role, and hence, we expect 
S2 to be a little less dominant within the spin-wave continuum. 

The next feature we discuss in this work is the behavior of S4 as a function of k for fixed 
uj. Figs. 2a, 3a and 4a show the behavior of S4 for 00/71 = 0.45,0.5 and 0.75 respectively. 
Note that we have scaled S4 to appropriate units. Figs. 2b, 3b and 4b show the behavior of 
62 as a function of k for the same values of uj. Let us for example discuss the case uj/tt = 1/2 . 
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We see from fig. 3b that the function S 2 vanishes for (roughly) k/n < 0.8. Looking back into 
fig. 1, this corresponds indeed to the region outside the spin-wave continum for u/n = 1/2, 
i.e., k/n < 5/6. The function S 2 starts at k/ir = 5/6 with a large value and goes to a 
minimum at k/n = 1 . 

Fig. 3a shows that the function £4 has a somewhat similar behavior. From the figure, we 
read that S4 too is not vanishing for (roughly) 0.8 < k/n < 1 . From fig. 1, that corresponds 
within the ^-region to 0.84 < k/n < 1. But fig. 1 shows also that for u/n = 0.5, S4 may be 
non- vanishing in the interval < k/n < 0.16. That contribution doesn't appear on fig. 3a, 
presumably because S4 there is negligeable. This is confirmed in the case u/n = 0.45 which 
is close to the case uo/n = 0.5: there we see S4 having a very small contribution in the 
corresponding interval. From fig. 3a, we see that S4, like S 2 , starts from a large value at its 
lower boundary k/n = 0.84 and decreases while moving to k/n = 1. 

The case u/n = 0.75 is practically the same. Fig. 4b shows S 2 starting from zero at 
roughly k/n = 0.25, getting to a maximum and sharply dropping to zero a little further. 
It then starts sharply again from a large value a little after k/n = 0.7 and decreases to 
a minimum at k/n = 1. This is also consistent with fig. 1: the function S 2 starts to be 
nonvanishing for u/n = 0.75 at k/n = 0.2447. It stays nonvanishing until k/n reaches the 
value 0.2699. It remains identically zero until k/n reaches the value 0.7301 at which we enter 
back into the ^-region. 

As we said, S4 has in this case too the same overall behavior as that of S 2 . In fig. 4a, we 
see that S4 starts to increase from the value zero at k — 0. It goes quickly to a maximum 
before dropping sharply to zero a little before k/n = 0.3. It stays at zero till a little after 
k/n = 0.7 and increases sharply. Then it decreases while wiggling to k/n = 1. This overall 
behavior is also consistent with fig. 1. Indeed, for uo/n = 0.75, the ^-region starts at k — 
and extends first to k/n = 0.2413. Then we get outside this region from this value of k 
till k/n = 0.7587. Then S4 is no more identically zero beyond this point until we reach the 
point k = n. As we said, this is quite consistent with fig. 4a. 

V. DISCUSSION AND CONCLUSION 

In this work, we have discussed the behavior of the exact four-spinon contribution S4 
to the dynamical correlation function S of the s = 1/2 antiferromagnetic Heisenberg model 
and compared it to the one of the exact two-spinon contribution S 2 . We first reviewed the 
model and the spinon structure of the corresponding Hilbert space. We then gave a brief 
account of the results concerning S 2 and its comparison to the Miiller ansatz. The first thing 
concerning £4 we discussed is the region in the (k, u>)-plane in which £4 is not identically 
zero. We found it to be different and larger than that of S 2 , i.e., the spin-wave continuum. 
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Both regions are drawn in fig. 1. We then discussed the behavior of S4 as a function of k for 
fixed values of u, and compared it to the one of S2. These behaviors are plotted in figs. 2,3 
and 4. We have found that the overall shape of S4 is more or less the same as that of S2. 

The first thing to emphasize we think is the overall similarity between the shapes of S4 
and 5*2. This is not at all expected from the outset, given the more complicated expression 
of S4, see (P3"D, as compared to that of S2, see fllBp. Does this mean that the shape of the 
other n > 2-spinon contributions and hence of the total S is already more or less "traced" 
by that of S2? We think that at this stage, it is too early for such an inference: this is a 
preliminary investigation into the behavior of S4 and clearly more work is needed. 

In any case, it would have been interesting to measure S4 for other (larger) values of 
u, especially in regions where S2 is identically zero whereas S4 is not, see fig. 1. But 
as us increases, the structure of the function F(pi, ...,^4) of eq(P^) in the (p 3 ,p 4 )-plane gets 
"richer" , which means numerically harder to handle. To illustrate this point, we have plotted 
for the reader in figs. 5 the function F for k/n = 1/2 and 00 = 2tt (5a), uj = 3tt (5b). One 
can see that as u increases, the function F gets distributed nontrivially in larger areas in the 
(p3,P4)-plane, with an increasing more involved structure. This is the main reason why we 
preferred to defer the discussion of these regions to future work. Also, for the same reason, 
we have deferred the systematic discussion of S4 as a function of u for fixed values k. 

One other interesting question we haven't touched on in this work but merely alluded 
to at the end of section 3 is the following: how much 5*4 accounts for in the total S? In 
other words, is S2 + S4 better an approximation to the total S than S2 alone, and if yes, 
by how much? As we said, we know that 5*2 accounts for about 70% of the total S, which 
means that roughly 30% are left for the n > 2-spinon contributions. To tackle this question 
as it should, one has to rely on a certain number of sum rules S is known to satisfy exactly. 
Then one compares the contribution of S2 + S4 to the exact result and carries a discussion 
thereon. 

The other interesting question one may ask is the physical interpretation and implications 
of the behavior of the four-spinon DCF we have described in this work. What would certainly 
be interesting is to be able to systematically measure S outside the spin-wave continuum so 
that we are assured of having the effects of the two-spinon contribution eliminated. 
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Figure captions 



Fig. 1: The regions in the (k, c<j)-plane inside which S2 (dashed lines) and S4 (solid lines) 
are not identically zero. The boundaries are indicated as defined in the text. 

Figs. 2-4: S4 (a) and S2 (b) as functions of k for fixed to. The values of u are indicated. 
Note that for the sake of illustration, £4 is plotted to appropriate units. 

Figs. 5: The function F(pi, P2, P3, P4) of relation (|23|) plotted in the (^3, ^4)-plane for k = 
0.5tt and u = 2tt (a), uo = 3tt (b). Note that each F is scaled to appropriate units. 
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